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ABSTRACT 

We report exact results for one-dimensional reaction-diffusion models A + A —>■ 
inert, A + A — > A, and A + B — > inert, where in the latter case like particles coagulate 
on encounters and move as clusters. Our study emphasized anisotropy of hopping 
rates; no changes in universal properties were found, due to anisotropy, in all three 
reactions. The method of solution employed mapping onto a model of coagulating 
positive integer charges. The dynamical rules were synchronous, cellular-automaton 
type. All the asymptotic large-time results for particle densities were consistent, in the 
framework of universality, with other model results with different dynamical rules, when 
available in the literature. 
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I. INTRODUCTION 



Diffusion-Limited Reactions (DLR) involving aggregation and annihilation pro- 
cesses are important in many physical, chemical and biological phenomena [1-3] such 
as star formation, polymerization, recombination of charge carriers in semiconductors, 
soliton and antisoliton annihilation, biologically competing species, etc. In this paper 
we generalize and apply a recently introduced method [4-5], in order to study by exact 
solution effects of anisotropy in some common DLR in one dimension (1-D), specifically, 
A + A — > A, A + A — > inert, and a two-species annihilation model A + B — > inert in 
which like particles coagulate irreversibly. The detailed definitions will be given later. 

Scaling approaches and other methods have yielded [1-3,6-12] the upper critical 
dimension, D c , for various reactions. Typical values range from 2 to 4. For spatial 
dimensions lower than D c the kinetics of these reactions is fluctuation-dominated, and 
we cannot expect the rate equation approach to be valid. Indeed the mean-field rate 
equation approximation ignores effects of inhomogeneous fluctuations. Fluctuation- 
dominated DLR have been subject to numerous studies by other methods, notably, 
exact solutions and asymptotic arguments [1,13-22] in ID. The ID reactions have also 
found some experimental applications [23-24]. These studies have assumed isotropic 
hopping (reactant particle diffusion). 

Recently Janowsky [25] concluded, based on numerical results and phenomenologi- 
cal considerations for the two species annihilation reaction A + B — > inert, that making 
the hopping fully directed would change the universality class in ID. Specifically, the 
large-time particle concentration (assuming equal densities of both species), would scale 
according to c(t) ~ t -1 / 3 instead of the isotropic-hopping power law t -1 / 4 . A few exact 
and numerical results available in the literature on anisotropic reactions involving only 
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one species [26-27] indicate that the power law is not changed. The model of [25] as- 
sumed that like particles interact via hard-core repulsion; this seems to be an essential 
ingredient for observing the changeover in the universality class. 

In this work we report the exact solution for two-particle annihilation with anisotropic 
hopping. We consider discrete-time simultaneous-updating dynamics, also termed syn- 
chronous dynamics, so that our models are cellular-automaton type. The universality 
classes of behavior at large times and large spatial scales are expected not to depend 
on details of the model dynamics. However, in order to achieve exact solvability we 
took "sticky-particle" rather than hard-core interactions: the like particles coagulate 
on encounters and diffuse as groups. Our exact calculations yield the t -1 / 4 power law, 
found earlier for similar "sticky-particle" reactions with isotropic hopping [28-29] . This 
result is probably due to absence of hard-core interactions in our model. 

For unequal initial concentrations, the large-time behavior changes, see [28-29]. The 
crossover between the two regimes is derived analytically. Finally, we also obtain new 
exact results for single-species two-particle aggregation and annihilation reactions with 
anisotropy; see also [26] . We find that anisotropy does not change the universality class 
of kinetics of these reactions. A short version of this work was reported in a letter-style 
publication [30]. 

This paper is organized as follows: in Section II we define and review some existing 
results for the models studied. In Section III our method for exact solution is introduced. 
In section IV, results for the single-species models are presented. Finally, Section V is 
devoted to the two-species model. Further scaling analysis of the two-species model and 
the brief summarizing discussion are left to Section VI. 
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II. REACTION-DIFFUSION MODELS 



In lattice DLR models with like particles it is usually assumed that the particles 
hop independently, to the extent allowed by their interactions, to their nearest neighbor 
sites. Whenever two particles meet, they can both annihilate which corresponds to the 
reaction A + A — > inert. If, however, only one particle of a pair disappears we get the 
reaction A + A — > A usually termed "aggregation." Such single-species reactions have 
the upper critical dimension D c = 2 [6-9] . For D < 2 the particle concentration at large 
times behaves according to c(t) ~ t~ D l 2 . Specifically, the ID kinetics of these reactions 
is no n- mean-field, with the typical diffusional behavior c(t) ~ t -1 / 2 . This result is not 
affected by short-range interactions between the particles and is not sensitive to their 
initial distribution as long as initial correlations are sufficiently short-range. 

Consider now the two-particle annihilation model, to be termed the AB model. 
Particles hop randomly to one of their nearest-neighbor sites. Whenever two parti- 
cles meet, unlike species annihilate, A + B — > inert. When like species meet, some 
interaction must be assumed. The simplest interaction is hard-core: diffusion attempts 
leading to multiple occupancy of lattice sites are discarded. For such reactions, assum- 
ing equal average A- and .^-concentrations and uniform initial conditions, the upper 
critical dimension is D c = 4 [6-9]. The (equal) particle concentrations scale according 
to c(t) ~ t~ D ^. A surprising largely numerically-based recent result of Janowsky [25] 
is the new exponent ~ 1/3, replacing 1/4, for anisotropic particle hopping in ID. For 
unequal initial concentrations, the density of the minority species is expected to decay 
faster than the symmetric-case power-law; some specific results will be referred to later. 

In order to obtain a solvable model in ID, let us now consider the AB annihilation 
model with the "sticky particle" interaction. Thus, like particles coagulate irreversibly 
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on encounters, for instance, 

nA + mA — > (n + m)A , (2.1) 

and the clusters thus formed then diffuse as single entities with the single-particle hop- 
ping rates. When unlike clusters meet at a lattice site, the outcome of the reaction 
is 

{(n — m)A if n > m 
inert if n = m (2.2) 

(m — n)B if n < m 

Recent numerical results and scaling considerations for these reactions [31] in D = 1, 2, 3 
seem to suggest that they are mean-field in D = 2, 3. However, in ID the power-law 
exponent for the density is 1/4, obtained by exact solution [28-29] which also yielded a 
faster power-law decay ~ t -3 / 2 for the minority species in case of unequal densities of 
A and B. 

As mentioned earlier, in order to obtain exact solutions for our single-species and 
AB models, we first solve another model [4], of coagulating nonnegative integer charges 
in ID [4-5], with anisotropic hopping. All our models are defined with synchronous 
dynamics to be described in detail in the next section, along with detailed dynamical 
rules and their exact solution. From the coagulating-charge model one can derive results 
for the reactions A+A — > A and A+A — > inert with anisotropic hopping; see Section IV. 
Some of our expressions are new, while other are consistent with results available in the 
literature. We next use the approach of [28-29] to extend the coagulating-charge solution 
to the "sticky" AB model. We find that hopping anisotropy does not change the density 
exponent in ID: it remains 1/4; see Section V. 
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III. SYNCHRONOUS DYNAMICS OF COAGULATING CHARGES 



Let us consider a one-dimensional lattice with unit spacing. Following [4], we 
consider diffusion of nonnegative charges on this lattice. Initially, at t = 0, we place 
positive unit charge at each site with probability p or zero charge with probability 

1 — p. Furthermore, we consider synchronous dynamics, i.e., charges at all lattice sites 
hop simultaneously in each time step t — > t + 1, where t = 0, 1, 2, . . .. However, the 
probabilities of hopping to the right, r, and to the left, £ = 1 — r, are not necessarily 
equal. Since this dynamics decouples the even-odd and odd-even space-time sublattices, 
it suffices to consider only those charges which are at the even sublattice at t = 0. Thus 
we only consider the lattice sites j = 0, ±2, ±4, ... at times t = 0, 2, 4, . . ., and lattice 
sites j = ±1, ±3, ±5, ... at times t = 1, 3, 5, . . .. 

One can view the diffusional hopping as taking place on a directed square lattice 
with the time direction along the "directed" diagonal. This lattice is illustrated in 
Figure 1; note that the charges arriving at site j at time t can only come from the sites 
j — 1 and j + 1 at time t—1. Finally, the "interaction" between the charges is defined by 
the rule that all charge accumulated at site j at time t coagulates. There can be 0, 1 or 

2 such charges arriving from the two nearest neighbors of j in the time step t — 1 — > t, 
depending on the random decisions regarding the directions of hopping from sites j ± 1 
in this time step. 

This model can also be viewed as diffusion-coagulation of unit-charge "particles" 
C, where the coagulation is represented by the reaction 



nC + mC ^ (n + m)C . (3.1) 

Such reactions, without the limitation of positive or integer charges, and with an added 
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process of feeding-in charge at each time step, with values drawn randomly at each site 
from some fixed distribution, have been considered as models of self-organized criticality 
and coagulation [5,32-33]. These studies were limited to isotropic hopping, i.e., r = I = 
i. Our interest in these reactions is in that their dynamics can be mapped [4,28-29] 
onto that of both the single-species (Section IV) and "sticky" AB (Section V) reaction- 
diffusion models introduced in Section II. However, before discussing and utilizing this 
mapping, let us present the exact solution of the model of coagulating charges with 
anisotropic hopping, following the ideas of [4-5]. 

For each time t and at each lattice site j (of the relevant sublattice) we define 
stochastic variables, 



- (+\- ! Probability r f , 

Tj[T) ~\0, probability £ {6 } 

which represent the hopping direction decisions. Then the stochastic equation of motion 
for the charges qj(t), equal to the number of C particles at site j at time £, is 



q n (t+l) = T n - 1 (t)q n - 1 (t) + [l-T n+1 (t)]q n+1 (t) . (3.3) 

The total number of C-particles, or the total charge, in an interval of k consecutive 
proper-parity-sublattice sites, starting at site j at time t, is given by 



fc-i 

Sk,j(t) = ^q j+ 2i(t) = qj(t) + q j+2 (t) + ■■■ + +q l+2k - 2 (t) . (3.4) 

Due to conservation of charge in this process, the equations of motion (3.3) yield the 
following relation, 
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SkAt + !) = T n -i{t)q n -i(t) 

+ q n +l(t) H h q n+2 k-3(t) 

+ [l-T n+ 2k-i(t)]q n +2k-i(t) • (3-5) 

This result indicates that only the two random decisions at the end points are in- 
volved in the dynamics of charges in consecutive-site intervals. The exact solvability 
of coagulating-charge models is based on this property, as first observed in [5]. Other 
solution methods were used in related interface-growth models [34] which will not be 
discussed here. 

Let us introduce the function 

7(s,m) = S Sj7n 

and averages, 

f k , m (t) = {I(S k , j (t),m)) . (3.7) 

The average (• • •) is over the stochastic dynamics, i.e., over random choices of the 
decision variables r^t), as well as over the random initial conditions. Since the latter 
are uniform, the averages fk,m(t) in (3.7) do not depend on the lattice site j. Functions 
other than the Kronecker delta have been used for I(s,m) in the literature [4-5,32-33]. 
With our choice (3.6), the resulting averages fk,m(t) correspond to the probability to 
find exactly m charge units in an interval of k consecutive sites. For instance, fi, m (t) 
is the density (fraction) of sites with charge m. 

Note that (3.5) essentially represents the following simple rule: charge "fed" into 
an interval of k sites comes from k + 1 sites, at the preceding time-variable value, with 
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(3.6) 



probability rl, from k — 1 sites with probability rZ, and from two possible groups of 
k sites, with probabilities r 2 and £ 2 ; this is illustrated in Figure 2. Furthermore, the 
variables Ti(t) and Sk, n (t) are statistically independent because the latter only depends 
on "decision making" variables r« at earlier times. Therefore, the averages introduced 
in (3.7) satisfy, for any function J(s,m), the following equations of motion, 



fk,m(t + l)=r£ [fk+l,m(t) + fk-l,m(t)] + + /*,»»(*) • ( 3 - 8 ) 

Interestingly, the m-dependence is parametric in (3.8). However, it does enter the 
initial conditions. It is also convenient [4] to define fo,m(t) = 7(0, m) in order to extend 
the applicability of (3.8) to all t > 0. For our specific choices, we have the following 
expressions. Firstly, the initial conditions of placing charges 1 or at each site, with 
respective probabilities p and 1 — p, correspond to 



/^(0) = {f (1 - rt '"™^ °^.- k 

^ m > k 

where ( ) = k\j (m\(k — m)\). Secondly, the boundary condition is that a null interval 
cannot have charges, i.e., 



f ,m(t > 0) = 5 0) m • (3.10) 

In order to solve the equations of motion (3.8) we introduce the double generating 
function of fk,m(t), over the time variable t and over the number of charges m, with 
fixed number of sites k, 



oo oo 



g k (u, w) = J2J2 fk,m(tW«> m ■ (3- 11 ) 

t=0 m=0 
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It is also convenient to introduce the variable a = r — £ directly measuring the hopping 
anisotropy, 



r=(l + a)/2 and £=(l-a)/2. (3.12) 

A straightforward but tedious calculation then yields the following difference equation 
which derives from the equations of motion (3.8), with(3.9), 



(l + a 2 ) u - 2 4 k 

g k+1 (u,w) + 2— r 2n gk{u,w)+g k -i(u,w) = ^—(wp + l-p) . (3.13) 

(1 — a z ) u (1 — ar) u 

The initial and boundary conditions "translate" as follows, 



g k (0,w) = (wp + l-p) k , (3.14) 



g (u,w) = — *— . (3.15) 
1 — u 

The solution of (3.13) is obtained as a linear combination of the special solution 
Q(wp+1— p) k proportional to the right-hand side, and that solution of the homogeneous 
equation which is regular at u = 0. The coefficient Q is obtained by substitution, 



O = 4(wp + l-p) 

(1 - a 2 )u(wp + 1 -p - A+) (wp + 1 -p- A_) ' 

where it is convenient to express the denominator in terms of the roots A± of the 
characteristic equation of (3.13), 

, 9 (l + a 2 ) u - 2 , 

A 2 + 2 A- ^- A + 1 = 0. 3.17 

(1 — a 1 ) u 



These roots are given by 



2- (l + a 2 )u±2^/(l — u) (1 — a 2 u) 



(3.18) 



(l-a 2 )u 



and the root A_, which is nonsingular as u — > 0, also gives the homogeneous solution 
proportional to . The proportionality constant is determined by (3.15). In summary, 
the solution takes the form 



The dependence of this result is via Q, see (3.16), and it is of a simple rational 
form. Therefore, derivation of the dependence on the number of charges, "generated" 
by w, is relatively simple. In our applications we will concentrate on densities of re- 
actants at single lattice sites, derived from fk=i,m(t)- The m-dependence then follows 
by expanding (3.19) in powers of w. However, the -u-dependence of (3.19) with (3.16), 
(3.18), is more complicated. Therefore we will use the generating functions for the 
time-dependence and most of our explicit time-dependent expressions will be derived 
as asymptotic results valid for large times. Indeed, the power series in u are then con- 
trolled by the singularity at u = 1, and analytical results can be derived by appropriate 
expansions. Specifically, let us introduce the time-generating function for the quantities 
fi,m(t) which represent the probability to find charge m at a lattice site at time t. We 
define 




(3.19) 



oo 




(3.20) 



t=o 



The central result of this section is thus 
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G m (u) = 5 m fi 



1 — u (1 — a 2 ) u 



4A + p T 



(1 -a 2 )u(l -p- A+) 



ra+l 



(3.21) 



where A± were given in (3.18). Note that G m (u) is just the m th Taylor series coefficient 
in w, of the function gi(u,w). 
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IV. SINGLE-SPECIES REACTIONS 



In this section we map the coagulating-charge model onto the two single-species 
reactions introduced in Section II. This approach follows recent work [4], although re- 
lated ideas have been used in earlier literature, for instance, in [35]. Let us consider 
first the two-particle aggregation reaction A + A — > A. In the coagulating-charge model 
we now regard each "charged" site as occupied by an A-particle, and each "uncharged" 
site as empty of A-particles. Specifically, any charge m > represents an A particle. 
No charge at a site, m = 0, corresponds to absence of an A particle. The dynamics of 
the coagulating charges then maps onto the dynamics of the reaction A + A — > A on the 
same lattice, with the initial conditions of placing particles A randomly with probability 
p, or leaving the lattice sites vacant with probability 1 — p, so that the initial density of 
the A particles is 



We now observe that the quantity fi,o(t) gives the density of empty sites in both 
models. Therefore, the particle density (per lattice site), c(t), in the aggregation model, 
is given by 



c(0)=p . 



(4.1) 



c(f) = 1 - /i >0 (f) • 



(4.2) 



The generating function is therefore easily derivable from (3.21), 



oc 



1 



1 - A 



4(1 -p) 



E(«) = X)c(ty 



1 — u 



G (u) 



1-u 



+ 



(l-a 2 )u(l-p- A+) ' 



(4.3) 



t=o 
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The function E(u) is actually regular at u = 0. Thus the Taylor series is controlled 
by the singularity at u = 1, near which we have the leading-order term 



E(u) = 



vT 



u 



+0(1) 



(4.4) 



This yields the leading-order large-time behavior, 



C(£) « ; 

While we are not aware of other exact solutions for this model with anisotropic hopping, 
the result (4.5) is not surprising. Indeed, the leading-order large time behavior is ex- 
pected to be universal in that it does not depend on the initial density p. Furthermore, 
the diffusion constant T>(a) = (l — a 2 ) V(0) decreases proportional to 1 — a 2 when the 
anisotropy is introduced (single-particle diffusion is then of course also accompanied by 
drift). Therefore, as a function of V(a)t, the result (4.5) also does not depend on the 
anisotropy and in fact it is the same as expressions found for other A + A — > A reaction 
models, with different detailed dynamical rules (with isotropic hopping), e.g., [20]. 

Let us now turn to the two-particle annihilation model A + A — > inert. The appro- 
priate mapping here is to identify odd charges with particles A and even charges with 
empty sites. Indeed, the dynamics of the coagulating-charge model is then mapped onto 
the reaction A + A — > inert. Thus, each lattice site with an odd charge q = 1, 3, 5, . . . 
will be replaced by a site with one A particle. Each lattice site with an even charge 
q = 0, 2, 4, . . . is empty of A-particles. The rules of charge coagulation then reduce 
to the desired reaction. Furthermore, relation (4.1) applies. However, the generating 
function E(u) for particle density is given by a different expression, 
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(4.5) 



oo . . 

E(u) = G 2j+ i (u) = = ±? j . (4.6) 



M> 3 (l-ai)u[(l-p-A + ) 2 -p2] 



The large-time behavior is similar to the aggregation reaction, with the universal ex- 
pression which only differs from (4.5) by a factor of 2, 

c(t) K ■ (47) 

While the large-time behavior of both models is model-independent and otherwise 
universal as described earlier, the finite-time results do depend on details of the dynam- 
ical rules. For our particular choice of synchronous dynamics on alternating sublattices 
(Figure 1), there exists an exact mapping relating the isotropic-hopping aggregation 
and annihilation reactions [4] . This mapping was also found for the anisotropic-hopping 
results obtained here. Specifically, we find (by comparing their generating functions) 



2c ineTt (t;p) = c A (t;2p) , (4.8) 

where the subscripts denote the outcome of the reaction while the added argument 
stands for the initial density. Thus if we consider an annihilation reaction with, initially, 
at t = 0, half the particle density as compared to an aggregation reaction with the same 
synchronous dynamical rules, then the density ratio will remain exactly 1/2 for all later 
times, t = 1, 2, . . ., as well. 

We note that due to complexity of expressions involved, results like (4.4) and es- 
pecially (4.8) in this section as well as many other expressions in the next two sections, 
could only be derived by using symbolic computer programs Maple and Mathematica. 
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V. TWO-SPECIES ANNIHILATION MODEL 



In this section we present results for the AB model defined in Section II. We assume 
that initially particles are placed with density p, but now a fraction a of them are type 
A, and a fraction f3 are type B. Clearly, 



Furthermore, the initial A- and B-particle concentrations are, respectively, ap and (3p. 
The concentration difference is constant during the reaction; it remains (a — (3)p. At 
large times, this is also the limiting value of the density of the majority species, while 
the density of the minority species vanishes. In what follows we assume 



which can be done without loss of generality. Indeed, the results for a < f3 can be 
obtained by relabeling the particle species. Thus, either the concentrations are equal or 
the majority species is always A. Our goal will be to calculate the density, c(t), of the 
majority species A. 

The dynamics of the AB model can be related to that of the coagulating-charge 
model of Section III by adapting the ideas of [28-29]. First, we note that the dynamics 
of the "sticky" A + B — > inert model can be viewed as coagulation. Thus, if a group of n 
particles A and m particles B meet, they can be viewed as coagulating and continuing 
to move together. Of course, at later times this combined group may become part of 
a larger coagulated cluster of particles. However, if their contribution to the particle 
count is according to (2.2), the annihilation is properly accounted for by this counting 



a + (3=l . 



(5.1) 



a>0, 



(5.2) 



both upon original and later coagulation events with other clusters. Alternatively, one 
can view these particles as new charges, +1 for A, and —1 for B. If the total charge 
of a coagulated cluster is positive than we view it as a group of A particles (equal in 
their number to the charge value). If the charge is negative, we consider the cluster 
B-particle, while if the charge is 0, we regard this cluster as nonexistent (inert) in the 
AB model. 

The probability of having an m-particle (charge m) cluster in the original positive- 
charge-only model was given by /i, m (t). Each such m-particle cluster can have charge 
n = —m, — m+2, . . . , m— 2, m, where we now refer to the new, ± charge definition rather 
than to the positive charges of the original coagulation model. The key observation 
is that having a "species" label assigned to a particle at time t = is statistically 
independent of its motion and affiliation as part of clusters at later times. Statistical 
averaging over the particle placement initially, and their motion at later times while 
coagulating to form particle clusters, is uncorrelated with the choice, with probabilities 
a and /?, of which species label to allocate to each particle at time t = 0. 

Thus the density (per site) of m-size clusters with exactly n units of charge, where 
m-size is that of the original coagulation models, while the charge — m < n < m is the 
new ± type, is given by 

m-\-n m — n Tfl\ 

*mAt) = « 2 P 2 (m+ny (r»-n), A,™(*) ' ( 5 ' 3 ) 

Therefore the density per site of A-particles, i.e., the density per site of the + charge, 
can be written as 




(5.4) 
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Similar to Section IV, let us denote the time-generation function of c(t) by E(u); 
see (4.3). By using results from Section III, and after some algebra, we get the following 
expression for the generating function, 



E{U) = (l-a*)u( P + +K-l) " V h ] S(X ' V) • (5 - 5) 



Here we introduced the function 



s(x,y) = j2Y, xn+j y j [ n > ( 5 - 6 ) 

n=l j=o \ 3 J 



and the variables 



y = f • (5.8) 

To make progress, we have to evaluate the double-sum in (5.6). This is accomplished 
as follows. First, we write the sum 

oo 

S(x,y) = J2 s nX n , (5-9) 

n=l 

where 

_ " (xyY T(2j + n + 1) 

i!ro- + n + i) • (5 - 10) 



Now, by using the duplication formula for the gamma function, we have 



where = T(z + j)/T(z). The sum in s n is a special case of the hypergeometric 
function, 



s n = 2F 1 (v,v+ 1/2; 2z/;C) , (5.12) 

where z/ = (n + l)/2 and ( = Axy. 

Fortunately, this can be expressed in elementary terms; all subsequent references 
in this paragraph are to formulas in Chapter 15 of [36]. First, using Gauss' linear 
transformation (15.3.6), 

*n = T^^M " + V2; 3/2; 1 - 

+ r^rg+i/V 1 " 0" 1/2 » ~ 1/2; 1/2; 1 - C) • (5-13) 

Next, from the recursion relation (15.2.20), 

2 F 1 (u, v - 1/2; 1/2; 1 - C) = C 2 F 1 (u, v + 1/2; 1/2; 1 - C) 

+ (1 -2u)(l-C) 2 F 1 (v,v+l/2; 3/2; 1 - C) • (5.14) 



Now, by (15.1.10), 
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2 F! + 1/2; 3/2;l-C) 



(i-C)- 1/2 

2(1 - 2i/) 



(l + v^) 1 ^-(i-a/w) 



l-2z/ 



while (15.1.9) yields the analytical expression, 



(5.15) 



2 F 1 (i/,i/ + 1/2;1/2;1-C) = 2 



(i + v/w) 2 " + (l - >/i^c) 



-2f 



(5.16) 



and so we have 



2 F 1 (u, v - 1/2; 1/2; 1 - Q = \C [(l + V 7 !^) 2 ' + (l - V^X) 

+ l(i - 1/2 [(l + v/w) 1 " 2 ^ - (i - v/W) 



-2i/ 



l-2i/ 



(5.17) 



Putting all this together we find 



Sn = 2 2 "- 1 (i - C)" 1/2 (i + a/w) 1_2I/ = (i - ^vY 1 ' 2 ( 



1 + V 1 - 4xy 



, (5.18) 



and finally, 



5(x, ?/ ) = (l-4x ?/ )- 1 / 2 ^f T 

n=l ^ 



2x 



2a; 



+ VI - 4zy / V 1 - 4 ^2/ (1 - 2x + VI - 4xy) 

(5.19) 



It is useful to introduce the parameter 6 = a — (3 > which measures the excess of 
A over B at time t — 0, 
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a=(l + &)/2 and p=(l-b)/2 



(5.20) 



Consider first the equal concentration case 6 = 0. The large-time behavior of the 
concentration c(t) is governed by the singularity at u = 1 of the generating function 
E(u). The form of the latter was evaluated near u = 1 from the expressions derived in 
this section, with the result, 



E(u) = 



(1-w) 3 / 4 



1 — p 



2(l-a 2 ) 1/4 VHl-a 2 ) 374 



[1 - u) 1 ' 2 + 0(1 - u) 



. (5.21) 



The leading-order behavior of the A-particle concentration follows from the first term 
in (5.21), while the second term will be further discussed in Section VI. We get 



c(t) 



(5.22) 



2r(3/4)(l-a 2 ) 1/ V/4 

The most significant feature of this result is that, similar to the single-species reactions 
considered in Section IV, the anisotropy, a, dependence can be fully absorbed in the 
diffusion constant, in terms of V{a)t = (l — a 2 ) V(0)t. The exponent 1/4 was derived 
in [28-29] for different (isotropic) dynamical rules. 



A similar expansion for fixed b > yields 



bp 1-b 2 2 (1 - b 2 ) (2 - b 2 p) 



*(«> = Tfi + V^Wi - (i -..)'/>/ ^ + 0(1 " U> ' <5 ' 23) 

The leading term in (5.23) corresponds to the constant contribution c(t) = bp + . . . 
which is expected since A is the majority species. In fact, expansions near u = 1 are 
nonuniform in the limits b — > + and b — > 0~. In deriving (5.23) we used for the first 



time the fact that the majority species is A. The approach to the constant asymptotic 
density value is given by the third term in (5.23), 



1 - b 2 ) 2 - b 2 p) , 
c(*)-6p« * ^ ^ff- . 5.24 

Note that this difference is just the density of the minority species B. As before, the 
anisotropy dependence of this leading-order power-law correction is fully absorbed in 
the diffusion rate, while the exponent is consistent with the results of [28-29]. Details 
of the crossover in the limit 6^0 will be discussed in the next section. 
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VI. CROSSOVER SCALING IN THE TWO-SPECIES REACTION 



As emphasized in the preceding section, the limit u — > 1 is nonuniform at b — 0, i.e., 
the pattern of the asymptotic large-time behavior changes at equal A- and S-particle 
concentrations. It is of interest to explore this behavior in greater detail within the 
standard crossover scaling formulation. In this approach, one seeks a combination of 
powers of variables each of which vanishes in the limit of interest, such that this so-called 
scaling combination can be kept fixed in the double-limit. The appropriate choice is 
expected to yield a nontrivial variation of quantities of interest in the limit, as functions 
of the scaling combination. 

In our case, the appropriate scaling combination turns out to be proportional to 
b/ (1 — w) 1 / 4 , as determined by inspection of various limiting expressions. It proves con- 
venient to absorb certain constants into the precise definition of the scaling combination 

0", 



The time-generating function E(u) studied in Section V will be now analyzed in the 
double-limit 6^0 and u — > 1~, taken with fixed values of a. From expressions derived 
in Section V, one can obtain 



a = 



VP (1- a 2 ) 1/4 6/(1 -n) 1 / 4 . 



(6.1) 



E{u) « p 1 (l — a 




R(a) 



(6.2) 



where R is termed the scaling function. Note that the first two prefactors are constants 
in the limit of interest. However, the power b~ 3 is necessary to ensure scaling function 
values of order 1 for a of order 1. 



The scaling function R can be derived exactly, 



a 3 (a + a/4 + a 2 ) , 
R(a) = i ; . 6.3 

Note that it is analytic at a = 0, where a oc b. Thus, at the expense of introducing the 
nonanalytic factor in a which is power-law in 1 — u, see (6.1), we managed to "blow 
up" the regime of small b. The scaling limit provides, as usual, a better understanding 
of the crossover in the limit b — > 0. Note that for a <C 1 the following small-argument 
expansion of R(cr) applies, 

R(*) = \°* + \° A + O (a 5 ) . (6.4) 

It is interesting to note that the leading term here actually reproduces the first term in 
(5.21). The latter was the limiting form for u — > 1 at b = 0. Indeed, the 6-dependence 
cancels out, while the (1 — -^-dependence is the identical, simple power-law in both 
limits. However, the second term in (5.21) does not seem to correspond to the next 
scaling-expansion contribution; see (6.4). Corrections to the leading scaling behavior 
contribute to this term in the 6 = expansion. 

In the opposite limit, a — > +oo, we get the expansion 



R(a) = a 4 + 1 - 4a" 2 + O (a" 4 ) . (6.5) 

The first term here reproduces the leading term in (5.23). Indeed, the limit a — > +oo 
corresponds to u — > 1 at fixed small positive b. Interestingly enough, the next two 
terms in (5.23) are also reproduced in their small-6 form by the next two terms in (6.5). 
For instance, the second term in (6.5) yields 1/ [(l — a 2 ) b 3 p] in E(u). Similarly, the 
third term in (5.23) is reproduced with numerator 4 which is the correct small-6 limiting 



value. Thus the minority-species concentration (5.24) with, for small 6, the numerator 
replaced again by 4, is also contained in the scaling form. Of course, corrections to 
scaling, not discussed here, yield improved results. 

The main point of the scaling description is that it provides a uniform limiting ap- 
proximation in the double-limit 6^0 and u — > 1. Specifically, the region of nonuniform 
behavior near b = is exploded by the large factor ~ (1 — w) -1 / 4 . In terms of cr, the 
behavior is smooth and well defined. For instance, the result (6.3) applies equally well 
for a < which corresponds to A becoming the minority species. The limit of u — > 1~ 
at small fixed b < is described by the limit a — > — oo. The appropriate expansion 
takes the form 

R(a) = -1 + 4a" 2 + O (a" 4 ) , (6.6) 

similar in structure to (6.5) but without the constant-density first term. 

In summary, we derived exact results for several reaction-diffusion models in ID. 
The leading-order large-time particle densities show expected power-law and universal 
behaviors. Anisotropy of hopping has no effect on the universality class of the models 
studied, and it can be largely absorbed into the definition of the diffusion constant. 
While finite-time results are expected [26] to be more sensitive to the value of the 
anisotropy parameter a, they are cumbersome to derive and of less interest than the 
leading-order expressions. One interesting exception is the duality relation (4.8) which 
applies for all finite time values in our synchronous-dynamics models. 
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FIGURE CAPTIONS 

Figure 1: The ID even-odd sublattices represented as the two-dimensional space-time 
lattice directed along the time axis. The solid-line bonds show possible hopping 
event directions. 

Figure 2: The charge in a continuous span of k lattice sites at time t + 1 can come 
from (a) k — 1 sites at time t, with probability r£. Note that the directions of 
hopping from these k — 1 sites, shown by double-arrows, are immaterial. Only the 
two exterior sites of the larger, (k + l)-interval shown determine the probability 
r£. Another possibility is (b) for the charge to come from k sites at time t. In this 
case both end-sites of the (k + l)-interval hopped to the right. The probability of 
(b) is therefore r 2 . Similarly, the probability of the charge coming from the other 
fc-interval, event shown as (c), is £ 2 . Finally, (d) the charge can also come from all 
k + 1 sites shown at time t, with probability r£. 
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